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Abstract 

The stationary states of systems with many molecular motors are studied theo- 
retically for uniaxial and centered (aster-like) arrangements of cytoskeletal filaments 
using Monte Carlo simulations and a two-state model. Mutual exclusion of motors 
from binding sites of the filaments is taken into account. For small overall motor con- 
centration, the density profiles are exponential and algebraic in uniaxial and centered 
filament systems, respectively. For uniaxial systems, exclusion leads to the coexistence 
of regions of high and low densities of bound motors corresponding to motor traffic 
jams, which grow upon increasing the overall motor concentration. These jams are in- 
sensitive to the motor behavior at the end of the filament. In centered systems, traffic 
jams remain small and an increase in the motor concentration leads to a flattening 
of the profile, if the motors move inwards, and to the build-up of a concentration 
maximum in the center of the aster if motors move outwards. In addition to motors 
density patterns, we also determine the corresponding patterns of the motor current. 
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Introduction 



Cytoskeletal motors such as kinesin, dynein, and myosin are proteins which convert the 
chemical free energy released from the hydrolysis of adenosine triphosphate (ATP) into 
directed movements along filaments of the cytoskeleton. In cells, these motors drive var- 
ious transport processes, and are also involved in cell divisio n , cell locomoti on, and re- 



organization of the cytoskeleton (Sc hliwa and Woehlke 



2003 



Howard, 



20011 ). A lot of 



knowledge has been obtained from in vitro motility assays which allow for the measure- 
ment of single motor properties such as their y elocities, average walking distances, step 



2001). These quantities 



sizes, and the forces they exert (iHoward . 

for various types o f processive motors inc l uding conventional kine sin 



rave been measured 



Block et al 



1990; 



Svob gda et al 



Schnitzer and Block 



1993 



1993), Myosin V 



Mcvhofer and Howar< 



Mehta et al., 



1999 



cessive monomeric kinesin KIF1A (Okada and Hirokawa 



Howard ct al 



Veigel et al 



1999 



1989; 



1996; 



Vale et al 
2002). the pro 



Tomishige et al 



King and Schroer . 



2002), 



2000h . These motil- 



and cytoplasmic dynein (|Wang and Sheeta . 200ol : 
ity assays study systems consisting either of (i) mobile motors and immobilized filaments 
or of (ii) immobilized motors and mobile filaments. In addition, systems where (iii) both 
motors and fila ments are mobile and filaments can b e displaced by motor s have also be 



studied (sec e.g. 


Takiauchi 


1991 


2001; 


Kruse and Jiilicherl ^ 


!00(l). 



Urrutia et al 



1991; 



Nedelec et al 



motor 


3 have also be 


. 1997: 


Surrev et al.. 



In all of these systems, motors and filaments interact via hard core interactions arising 
from their mutual exclusion. Indeed, both motors and filaments occupy a certain spatial 
volume which cannot be occupied by another molecular structure. In particular, motors 
bound to filaments exclude other motors from the binding sites of the filaments. T he latter 



exlusion effects were first addressed in our previous work ([Lipowskv et al 



2001) in which 



we introduced a general class of driven lattice gas models for this purpose. 

In the following, we use these driven lattice gas models in order to explore how the 
arrangement of the filaments affects the motor transport in closed compartments. We 
consider uniaxial and centered filament arrangements and present results for the station- 
ary patterns of both motor density and motor current. Both types of arrangements are 
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accessible to in vitro experiments and mimic structures of the cytoskeleton as observed in 
vivo. The uniaxial systems mimic the geometry of axons or fungal hyphae, while centered 
systems are realized, for example, in the aster-like structures of microtubules extending 
from centrosomes. For the uniaxial systems, we have previously shown that traffic jams 



build up easily as a consequence of mutual exc 



usion (ILipowskv et al 



200J), while pre- 



2001) did not incorporate this mutual 



vious work on centered systems (N edelec et al 
exclusion. 

We will show in the following that uniaxial and centered systems exhibit rather differ- 
ent jamming behavior. While in uniaxial systems jammed regions grow upon increasing 
the motor concentration and spread over the whole system, the effect of jamming in cen- 
tered systems is less dramatic and jams remain small in this case. Increasing the motor 
concentration, however, influences the density profile in the non-jammed region. In addi- 
tion, we show that the traffic jams in uniaxial systems are rather insensitive to the motor 
behavior at the end of the filaments. In contrast, the latter behavior is crucial for the 
presence of jams in centered systems. 

The density profiles discussed here theoretically can be directly measured in biomimetic 



experiments in vitro, and, in fact, such de nsity profiles 



r ave re cently been measured for 



the case of centered or aster-like systems (|Nedelec et al 



2001). However, the latter ex- 



periment did not address the jamming behavior, which could be done by increasing the 
motor concentration in these systems. In addition, our theoretical density profiles can be 
compared to motor density profiles measured for the corresponding systems in vivo. Such 
in vivo d ensity profi 



systems. 



Seiler et al 



. e s hay e been reported for fungal hyphae, which represent uniaxial 
(2000) have observed motors localized at the tip of these hyphae, 
which corresponds again to the case of low motor density. In vivo, the motor concentra- 
tion can be changed by changing the level of expression of the corresponding gene; in that 
way ja m-like density p r ofiles have recently b een observed for another fungal kinesin-like 



motor (Konzack 



2004 



Konzack et al 



2005). The effect of exclusion (and, thus, jam- 



ming) is enhanced if the motors transport large cargoes such as membranous organelles. 
Jam-like behavior of organelles has been observed in axons (W. Saxton, private commu- 



3 



nication), extreme cases induced by mutations of motors (which are lethal in later 



of development) are accompanied by strong swelling of the axon (|Hurd and Saxton , 



Martin et al 



stages 



1996; 



1999). 



Our article is organized as follows. We introduce the theoretical model in the following 
section. In the sections 'Density profiles for uniaxial filament systems' and 'Density profiles 
for centered filament systems' we discuss jamming effects in two types of filament systems 
and present results for the motor density patterns obtained from Monte Carlo simulations 
and from a two-state model. Finally, we relate our results to recent experiments in the 
discussion section. The appendices describe the theoretical methods used in this article 
and some analytical calculations. 



Lattice models for molecular motors and filaments 



In this article, we study the stationary profiles of the motor density which build up within 
closed compartments containing filaments. These s tationary states are c haracterized by 



the balance of bound and unbound motor currents ([Lipowskv et al 



200ll ). Unbinding of 



motors from the filaments reflects the finite binding energy of the motor-filament com- 
plex which can be overcome by thermal fluctuations and leads to peculiar random walks 
of the motors which consist of alternating sequence s of directed motion along filaments 



and non-directed diffusion i n the surrounding fluid l|Aidari 



Nieuwenhuizen et al 



2002, 



1995 



Lipowskv et al 



2001; 



20041 ). see Fig. Ufa). In order to stud y these random walks, 



we have recently introduced lattice models ([Lipowskv et al 



200 lj). One useful feature of 



these models is that one can incor porate motor-m o tor interactions such as the mutual 



exclu sion in a rather natural way ((Lipowskv et al 



2001 



Klumpp and Lipowskv 



2003, 



2004). Motor-motor interactions are especially important on the filaments: motors are 
strongly attracted to filaments, so that the local density of motors on these filaments will 
typically be large even if the overall motor concentration is rather small. The importance 
of motor-motor interactions is further increased if motors accumulate in certain regions 
of closed compartments. 

Mutual exclusion of motors from binding sites of the filaments has two effects: (i) 



Binding of motors to the filament is reduced for those filament segments which are already 
occu pied by many motors. This e ffect is directly observ ed in decoration experiments (see, 



e.g. 



Song and Mande lkow , 



1993; 



Harrison et al 



1993). (ii) The mutual hindrance slows 



down the movement of motors in regions of high motor density. This second effect has 
not yet be en studied experim entally, but there are indications of it in microtubule glid- 



ing assays (Boh m et al 



2000). In addition, there is indirect evidence for such a slowing 



down from the self-organization of microtubules and motors, where an increase of mo- 
tor concentra t ion c an induce a transition from vortex to aster patterns of microtubules 



( Sur rey et al 



2001). From computer simulations, such a transition is expected if the mo- 



tors spend more time close to the end of a filament. This should happen if the motors 
are slowed down at the filament end by a traffic jam which builds up upon increasing the 
motor concentrations. 



Bound and unbound motor movements 



In the following, we describe the movements of molecu lar motors as random walks on a 



three-dimensional cubic lattice (Lipo wskv et al.1 



2001 



Nieuwenhuizen et al 



2005 



2004). 



One or several lines of lattice sites represent one or several filaments. The lattice constant 
I is taken to be the repeat distance of the filament which is 8 nm for kinesins moving 
along microtubules, so that filament sites of the lattice correspond to binding sites of the 
filament. A motor bound to a filament performs a biased random walk which describes 
the active movements along the filament. Per unit time r, it attempts to make forward 
and backward steps with probability a and (3, respectively. As backward steps are rare 
for cytoskeletal motors, we take = in the following which eliminates one parameter 
from our systems. Rather similar behavior is found for small nonzero values of (3. With 
probability 7, the bound motor makes no step, and with probability e/6, it unbinds to 
each of the four adjacent non-filament sites. The sum of all hopping probabilities per unit 
time r is one, i.e. the probabilities are related by 

a + /3 + 7 + 4e/6 = l. (1) 

When the motor particle reaches the end of the filament, it does not have the possibility 
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to step forward to another filament site. We will consider two different unbinding processes 
for this 'last' filament site as in our previous work (Klumpp and Lipowsky, 2003): 

- thermal unbinding: the motor particle detaches from the 'last' filament site with 
probability e/6 to the unbound site in the forward direction, but remains at the 'last' site 
with probability 7' = 7 + a — e/6, while the backward probability (3 and the sideward 
probability e/6 remain unchanged. Adjusting the no-step probability implies the modified 
normalization 

P + i + 5e/6 = 1 (2) 

for the hopping probabilities at the 'last' filament site. 

- active unbinding: the motor particle detaches from the 'last' site with probability 
a in the forward direction and with probability e/6 in the four sideward directions as for 
all other filament sites. In this case, the normalization of the hopping probabilities at the 
'last' filament site is given by (1). 

An unbound motor performs a symmetric random walk which corresponds to non- 
directed diffusive movement. It attempts to step to each adjacent lattice site with equal 
probability 1/6. If an unbound motor reaches a filament site, it can bind to this site with 
probability 7r a d. The random walk probabilities can be chosen in such a way that one recov- 
ers the measured transport properties of specific motors such as the bou nd state velocity, 



the unbound diffusion coefficient a nd the average walking distance (see 



2001; 



Klumpp and Lipowskv . 



Lipowskv et al 



2003). 1 The unbound diffusion coefficient D u ^ fixes the ba- 



sic time scale r = £ 2 /D u b. The probabilities a, 0, 7, and e are determined from the 

velocity = (a — j3)i/r of a single bound motor, the average walking distance along the 

filament Axb = 3^bT/(2e), the condition (3 = and Eq. ^ 

Mutual exclusion of motors is taken into account by rejecting all hopping attempts to 

lattice sites which are occupied by other motors. We take the motor particles to have a 
1 Note that the model used here does not account for the bound state diffusion coefficient or, equivalently, 
the randomness parameter of the motor movements. This para meter can be incorpora ted by introducing a 



second time scale for the movements of the bound motors fsee iLipowskv et all 1200 if) . Such an externded 
model leads to density profiles which are very similar to the ones described here. This indicates that the 
overall diffusion is essentially governed by the unbound diffusion process. 



6 



linear size comparable to the filament repeat distance I and to occupy a volume £ 3 . If the 
motors are attached to larger cargoes, exclusion is enhanced. In particular, a large cargo 
of linear size n£, when bound to the filament, effectively covers between nl and (2n — 1)1 
filament sites depending on the bound density. However, the functional relationships 



between the different densities and current are rather similar (MacDonald et al 



1968; 



1974). We will briefly discuss this case at the end of the paper 



McGhee and von Hip pel 
in the 'Discussion' section. 

These lattice models for systems with many molecular motors are related to driven 
lattice gas models which have been studied ext ensively in the context of non-equilibrium 



phase transitions (jKatz et al 



1983; 



Krug, 



1991). In the models studied here, the 'driving', 



i.e. the active directed movement, is restricted to the linear subspaces corresponding to 
the filaments. 

In the following, we will usually express all lengths and times in units of the filament 
repeat distance I and the basic time scale r, respectively. This means that the bound 
and unbound motor densities p\, and /9 u b that we will consider in the following are local 
particle number densities satisfying < /?b < 1 and < p u b < 1, which corresponds to 
< Pb < l/^ 3 and < p u b < l/^ 3 in dimensionful units. Dimensionful units will be used 
when presenting results for specific motor molecules. 



Filament arrangements and compartment geometries 

In this article, we study two types of filament arrangements within closed compartments 
as shown in Fig. ^b) and (c) . The first type is a uniaxial filament system where a closed 
cylindrical tube contains a number Nf of uniaxially arranged filaments, i.e. filaments ori- 
ented parallel to the cylinder axis and with the same orientation. We denote the coordinate 
parallel to the filament by x and the coordinates perpendicular to it by y and z. The tube 
has length L and radius R. 

The second type of system which we will study is a centered filament system, i.e., a 
radial or aster-like arrangement of filaments within a closed disk- like compartment. The 
number of filaments is again Nf. In this case, we denote the radial coordinate by r. The 
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linear extension of the compartment along the direction of the filaments, i.e. the disk 
radius, is denoted by L and the disk height by h. 

In both cases, we take all filaments to have the same length which is equal to the 
corresponding linear extension of the compartment, i.e. to the tube length and to the disk 
radius in the case of uniaxial and centered filament systems, respectively. Shorter filaments 
lead to very similar results. An example would be filaments in radial arrangements which 
are nucleated from a centrosome and extend from r = R c , the centrosome radius, to the 
disk radius r = R. In addition, since the compartments are closed, and the number of 
motors, denoted by N, stays constant within each compartment. 



Density profiles for uniaxial filament systems 



We first consider uniaxial arrangements of filaments within a closed tube as shown in 
Fig. mb). On the one hand, placing one or several filaments and motors inside a tube 
should be experimentally feasible. The tube could be either a glass tube as used for mi- 



cropipettes, a topographic channel as used for filament g uiding (ICle mmens ct al 



a liqu id microchannel on a chemically structured surface (|Gau et al 



1999; 



20031^ 



Brinkmann and Lipowskv 



2002). In all cases, tube diameters down to a few /^m can be achieved. On the other hand, 
tube-like geometries are also quite common in cells, the most prominent example being 
the axon of a nerve cell, a tubular cell compartment with a diameter in the range of 
few micrometers and a le ngth of up to a meter, which contains tens of microtubules per 



fim 2 ( Al berts et al 



2002J), i.e., typical distances of the microtubules are in the range of 
100 nm. Similar compartments, the hyphae, exist in the case of fungal cells. In addi- 
tion, some compartments inside the cell have tubular shapes and contain filaments such 
as strands of cytosol crossing vacuoles in plant cells, again with diameters in the micron 
range. 

We will now focus on the case of a single filament, since the case of Nf isopolar parallel 
filaments in a tube with cross-section <ft is essentially equivalent to a single filament in a 
tube with cross-section <f>/Nt provided that the filaments are equally distributed within 
the tube. (If the filaments are concentrated in a certain region, i.e., if the distance between 



filaments is small compared to the distance between filaments and the tube wall, depletion 
of motors is enhanced; depletion effects are rather weak, however.) Let us consider a 
cylindrical tube of length L and radius R with one filament located along its symmetry 
axis. Imagine now that a certain number of motors is placed into this tube. In the absence 
of ATP, the system attains an equilibrium state, where binding to and unbinding from the 
filament balance each other locally, i.e. at every single binding site. Both the bound and 
the unbound motor densities are constant and related by the radial equilibrium condition 

7T a dPub(l - Ph) = epb(l - Pub) ~ ePb, (3) 

where the terms (1 — pb) and (1 — p u b) describe mutual exclusion of bound and unbound 
motors, respectively, with (1 — p u b) ~ 1 for typical experimental situations. 

When ATP is added to the system, the motors start to move along the filament. We 
use the convention that the filaments are oriented in such a way that the bound motors 
move to the right. The motor current along the filament builds up a density gradient, 
which generates a diffusive current. In the stationary state, this diffusive current balances 
the drift current of bound motors. As a first approximation, we assume that Eq. El is also 
valid in the presence of ATP (which is justified if the velocity vy, is sufficiently small, as 
we will show below). The balance of currents can then be expressed by 

^bPb(l - Pb) = D uh <p ~ LW— ^-j^-. (4) 
dx vr ad dx 1 - p b 

It follows from this relation that, for low motor densities, the motor density increases 
exponentially along the filament and that motors accumulate at the right end of the 
filament, further increasing the importance of exclusion effects there. 

Simulation results 

Typical density profiles as obtained from Monte Carlo simulations are shown in Fig. El (a). 
If the total number N of motors is relatively small — one example is provided by A = 100 
in Fig. El — motors are essentially localized at the right end of the tube. Crowding of 
motors occurs only in a short region at the end of the filament where motors form a kind 
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of traffic jam. To the left of the traffic jam, the density has an exponential profile as 
predicted by the simple balance of active directed currents and diffusive currents. If a 
motor detaches from the filament in the crowded region, it will diffuse back over a certain 
distance and most likely rebind to the filament in the region to the left of the traffic jam. 
In the jammed region rebinding is strongly reduced, since essentially all binding sites are 
already occupied. Upon rebinding to the filament, the motor will move relatively fast to 
the right until it ends up in the jammed region again. 

These observations imply the coexistence of a low-density region with an exponential 
density profile and a crowded high-density region, separated by a relatively sharp domain 
wall or interface which corresponds to the beginning of the traffic jam. If the number 
of motors in the tube is increased, the jammed region spreads further to the left and the 
domain boundary is shifted towards smaller values of the spatial coordinate x, as shown in 
Fig. El for N = 350. Now, motors diffuse backwards over larger distances, since attempts 
to rebind to the filament fail, if the binding sites are already occupied. 

Finally, if the motor concentration is very large, there is only one domain with a high 
density of bound motors: The filament is crowded over its whole length and the bound 
density profile is essentially constant except for the regions close to the two ends of the 
filament, see the case N = 1000 in Fig. |2 In this case, motors may diffuse back over 
the whole system length, but both the diffusive current and the bound current along the 
filament are very small. 

The corresponding profiles of the bound motor current along the filament are shown 
in Fig. 12(b). The diffusive current of unbound motors integrated over the tube cross- 
section has the same absolute value, but the opposite sign. The current depends strongly 
on the position x along the filament as long as the filament is not completely jammed. 
Like the motor densities, it increases exponentially in the low density region. In the 
jammed region at the right end of the tube, the current decreases rapidly. It reaches its 
maximum close to the end of the traffic jam. Note that the maximal current accessible in 
these systems is smaller than t>b/4, i.e. smaller than the maximally possible current in a 
system with constant densities such as a tube system with periodic boundary conditions 
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(Kl umpp and LipowskvL l200ra ) . For the case of the completely jammed filament, the 
current profile is nearly flat, while the absolute current is small. 

To obtain a global characterization of transport in the system, we determined the 
average current defined by 

Jh = 7 / dxj h (x). (5) 
L Jo 

This quantity exhibits a maximum at an optimal motor concentration as a function of 
the total number N of motors in the system, i.e., as a function of the overall motor 
concentration, see Fig. |31 For small N, it grows linearly with the number of motors, while 
for large N it decreases again since motion of the bound motors is slowed down by the 
increasing traffic jam. However, this decrease of the currents is rather slow, since additional 
motors introduced into the system can only rarely find free binding sites. For the system 
shown in Fig. |21 the maximal current occurs for N ~ 350 motors, which corresponds to 
the intermediate case of the profiles in Fig. [21 

A second quantity, which gives a global characterization of the profiles, is the traffic 
jam length of the crowded domain. L* can be defined by the condition p\,(x*) = 1/2 
via = L — x*. Results for L* are also shown in Fig. El The three cases discussed above 
can now be distinguished as follows. For very small L*/L, crowding of motors only occurs 
in a small region at the filament end and the profile decays exponentially to the left over a 
large fraction of the system size. For intermediate values of L^/L with <C L*/L <C 1, the 
density profiles exhibit coexistence of domains with high and low bound motor densities. 
Finally, for L*/L ~ 1 the whole filament is crowded. Comparing the functional dependence 
of the traffic jam length L* with the average bound current Jt, shows that the optimal 
transport occurs when a large part of the filament is crowded, ~ 0.8L, but the traffic 
jam is not yet too dense. 

Two-state model 

To get some further insight into the properties of these self-organized density profiles, we 
studied the stationary states of these systems using a two-state model. In this model, 
which is described in detail in appendix the dependence of the concentration profiles 
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on the radial coordinate is neglected and motors can be in two states, namely bound and 
unbound. This approximation is justified, because the profile of the unbound motor density 
depends only weakly on the radial coordinate. The two-state approximation captures most 
of the relevant features of these profiles and numerical solutions for the stationary profiles 
are obtained much faster than by Monte Carlo simulations. 

For the uniaxial systems, the two-state model as defined in appendix^is given by 



vhPb{x) i - p h (x + 1; 



4>D 



ub 



Pub {x + 1) - Pub 0*0 



(6) 



and 



VbPb(x) 1-pbOc+l) -v h p h (x-l) l-pb(x) =Tf ad p uh (x) l-p h (x) -ep h (x) l-pub(x) ■ 

(7) 

These equations express the balance of bound and unbound currents and the binding to 
and unbinding from the filament, respectively, see appendix 

For the non-jammed low density region, some analytical results can be obtained from 
these equations which are presented in appendix[Bl In particular, an exponential increase 
of the density profile is obtained as pb ~ AAe X// with a length scale £ as given by Eq. IB.151 

In order to obtain results for arbitrary densities, we solved the two-state equations EH3 
numerically. Some profiles of the bound motor density as well as the average current as 
a function of the number of motors are shown in Fig. 0] for a relatively large system with 
L = 1000. While the main features are the same as for the smaller system discussed above, 
some additional details can be seen here. The current increases linearly with the number 
N of motors for small N, but at a certain point, N ~ 500 or N/L ~ 0.5 in Fig. E(b), 
the slope begins to change. The current then increases more slowly, but again nearly 
linearly, until it reaches its maximum. This change in slope of the current corresponds to 
the formation of a plateau in the density profile, where the density in the traffic jam is 
approximately constant and changes only little upon addition of motors. 

Until now, we have assumed thermal detachment, i.e., that detachment at the end of 
the filament occurs with the same rate as detachment at any other site of the filament. 
As mentioned before, see Eq. a second possibility is active unbinding, i.e., that mo- 
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tors detach with an increased rate at the end by making an active step which leads to 
unbinding. There is some indirect evidence for a quicker detachment of kinesin motors 
at the microtubule end from experiments and com puter simulation of the formation of 



aster and vortex patterns of microtubule by motors (S urrey et alJ . 120 01), In simulations, 
quick detachment at the microtubule end leads to the formation of vortex (or spiral-like) 
structures, while slow detachment is necessary for the formation of asters- like centered ar- 
rangements. Kinesin is able to form both asters and vortices, suggesting that detachment 



at the microtubule end is relatively quick, while the kinesin-relate d motor ncd 



o nly fo rms 



2001). 



asters and thus probably detaches slowly at the microtubule end (S urrey et al. 

We have determined density profiles for both cases using the two-state approach, see 
Fig. 0D These density profiles show that the jamming behavior is rather insensitive to 
the motor behavior at the end of the filaments. Except for the region very close to the 
filament end, the profiles for the two cases agree well. In particular, the domain wall 
or interface represented by the steep increase of the density profile at the beginning of 
the traffic jam is the same in both cases. This observation shows that the traffic jam is 
not due to the slow unbinding at the end, but due to the accumulation of motors with 
an exponential density profile, which follows from the balance of bound drift current and 
unbound diffusive currents in a uniaxial geometry. On the other hand, the density profile 
within the jammed region depends strongly on the detachment rate at the filament end. 
While there is a weak increase of the bound density inside the jammed region for small 
detachment rate at the filament end, the bound density decreases strongly in this region 
in the case of an increased detachment rate at the end, see Fig. |SJ 

Comparing the density profiles obtained from the two-state model with those obtained 
from simulations for the same parameter set, we find quite good agreement in the case 
where the entire filament is crowded by motors. For smaller overall motor concentrations, 
qualitative agreement is still good (except for the region close to the left boundary, see 
appendix but there are small quantitative discrepancies. We find that the length scale 
£ of the exponential increase of the density is smaller in the two-state model than in the 
simulations. Correspondingly the crowded region is slightly longer in this approximation. 
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This difference is due to neglecting a depletion zone close to the filament in the two-state 
model. Close to the filament, the unbound motor density is smaller than its radial average 
in the low-density region. Taking it as independent of the radial coordinate, we thus 
overestimate binding of motors to the filament. Therefore there are more motors bound 
to the filament in the two-state model than in the simulations, which results in a longer 
traffic jam, and the maximum of the current is shifted to a smaller number of motors in 
the system. 

In order to obtain a quantitative description of the radial density profiles, we have 
solved the full three-dimensional diffusion equation for the low-density region and derived 
an analytical expression for the depletion layer close to the filament as shown in appendix 
IU1 It follows from the latter expression that the radial profile of the unbound density is 
nearly flat far from the filament and exhibits a logarithmic depletion zone close to the 
filament. This confirms the observation that the unbound motor density depends only 
weakly on the radial coordinate which justifies the two-state approach. Comparing the 
results from this calculation with the simulation of the full model, good agreement is found. 
In Fig. El we have plotted both the longitudinal (a) and radial profiles (b) as obtained 
by both methods. The radial profiles exhibit the predicted depletion layer close to the 
filament in the low density region to the left of the traffic jam. In the crowded region, the 
unbound density is enhanced close to the filament in comparison to the value far from the 
filament. The full diffusion equation also leads to a condition for the length scale £ given 
by Eg. IC .291 and we obtain £ ~ 37.4 for the parameters used in Fig. in good agreement 
with the value from simulations which is £ ~ 37. In contrast, the two-state approximation 
yields the smaller value £ ~ 24, because it overestimates the current of motors binding to 
the filament. 

Density profiles for centered filament systems 

In this section, we consider profiles of the motor concentration in centered filament sys- 
tems or aster-like arrangements of filaments as shown in Fig. ^c). Such arrangements 
can be formed for microtubules in vitro either by nucleation from microtubule-organizing 
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centers tf Holv 



(Nedelec et al 



et all I1997T) or bv self-or ganization of microtubules and motor complexes 



1997; Su rrey et al 



2001). Centered filament systems mimic the most com- 



mon organization of microtubules in cells. Motivated by the restru cturing of this organi- 



zatio n during cell division and the formation of the mitotic spindle (jHvman and Kars enti 



1996), many experiments have focused on the case where the filaments are also mobile. 



In the following, we consider immobilized aster-like arrangements of filaments which 
are not reorganized by the action of motors. The asters consist of Nf filaments of length L 
arranged radially in a thin disk of radius L and height h. We take the filaments to extend 
from r = to r = L within the disk, but again smaller filaments lead to very similar 
results. In that case, active directed currents of motors along the filaments and diffusive 
motor currents will again be balanced in a stationary state. For the aster geometry some 
th eoretical and e x perim ental results for low motor densities have recently been reported 
by 



Nedelec et al 



( 200ll ). 2 We confirm their main theoretical result, an algebraic density 
profile far from the center of the aster, and extend the study of concentration profiles in 
asters by exploring the effect of mutual exclusion. 



Two-state equations for centered systems 

Centered filament systems are implemented in the two-state model by substituting the co- 
ordinate n used in the general expressions with the radial coordinate r, by using Eqs. IA.5T - 
IA.7l for the bound and unbound motor currents and introducing a geometrical weight factor 



In the experiments of iNedelec et al 



i 2001) the dynamics is more complicated, since they used mobile 
filaments and their motor constructs can also displace these filaments with respect to each other, so that 
the motor density profiles and the filament patterns develop in coordination. After some time, however, 
these systems reach a steady state, in which the filament pattern is stationary (although not completely 
immobile) and can, on average, be represented by a fixed filament system. In addition, once the aster-like 
structure is formed, the filaments are usually sufficiently separated from each other (with the exception of 
the center of the aster), so that the additional dynamics plays only a minor role. Indeed, in the low density 
case, the theoretical density profiles agree well with the experimental profiles, as we will discuss below. In 
order to exclude the more compicated dynamics, one could immobilize the filaments once the stationary 
state has been re ached or, alternatively, use centered microtubule systems nucleated from centrosomes 



( Holv et al 



1997) and conventional kinesins which cannot bind to two filaments at the same time. 
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4>(r) ~ for as described in appendix The latter factor implements the fact that with 
increasing r, the volume available for unbound diffusion increases. The two-state model 
equations are then given by 



VbPb(r) 1 - p b (r + 1) = -Dub <t>{r) Pub(r + 1) - Aib(j) 



(8) 



Vbpb(r) l-pb(r + l) -v h p h (r-l) l-p h (r) = 7fadPubW l-pb(r) ~ePb(r) l-Pub(r) 

(9) 

for the case of motors moving outwards in an aster. If motion of motors is directed inwards, 
i.e. if Vb < 0, the two-state equations are given by 



VbPb{r + 1) 1 - Pb(r) = 4>{r)D nh p uh (r + 1) - p uh (r) 



(10) 



Vbpb(r+1) l-pb{r) -v h pb(r) l-/9 b (r-l) = 7fad(r)p u b(r) l-pb(r) -ep h {r) 1-pubO") 

(11) 

In the lo w density limit , these equations lead to the algebraic density profile pb ~ r 71 



reported by 



Ncdelec et al 



(20011), see appendix El The exponent rj ~ Vb is positive if 



motors move outwards and negative if motors move inwards. 



Numerical results 

To study the effect of hard core exclusion in asters, we used the parameters given by 



Nedelec et al. for the motor constructs used in their experiments (Ne delec et al 



200LV ^ 



the numerics all parameters are given in units of the microtubule periodicity £ = 8 nm and 
the step time t s = 10 ms. Parameters of the bound state are: v b = 0.8/xm/s = £/t s and 
e = 0.01 r" 1 corresponding to unbinding after 100 steps, and those of the unbound state: 
D uh = 20/xm 2 /s = 3125 £ 2 /t s and 7^ = 2 .6 jmiWl 2 = 405.6 t ' 1 . Parameters which 



correspond to kinesin with beads as used by 



Lipowskv et al 



2001) lead to similar results. 



All results shown in the following are obtained for an aster of Nf = 300 microtubules of 
length 50 pm = 6250 £, which is confined in a slab of height 9/xm = 1125 £. 
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Motors moving inwards 

We consider first the c ase where motors mo ve inwards. For this case, experimental results 



have been reported by 



Nedelec et al 



(2001). Accumulation of motors in the center of the 
aster is observed by fluorescence microscopy. Profiles of the total motor concentration, i.e. 
the concentration averaged over bound and unbound motors, 

= ftWty>) i Pub „ i r ,-. + j^, (12) 

l + (j>(r) 4>{r) (p vr ad 

can be extracted from the fluorescence images. The last expression is valid for small motor 
densities and sufficiently large values of r and predicts that the density profile exhibits a 
crossover from a decay ~ r v_1 for small r to ~ r v for large r. This crossover behavior is 



seen in the experimental data ([Nedelec et al 



2001). 



For small overall motor concentrations, the numerical solution of the master equations 
exhibits the power law behavior predicted theoretically by neglecting exclusion effects. For 
the chosen parameters we find p u b{ r ) ~ with rj ~ —0.54 from the data for N = 10 
shown in Fig. a) in agreement with Eq. IB.23I In the center of the aster, a traffic jam is 
obtained already for small total number of motors. The traffic jam is however rather short 
and, in contrast to the case of uniaxially aligned filaments, does not grow substantially in 
length, when the number of motors in the system is increased, see Fig. Ufa). Jamming of 
motors occurs only for small r (< 20^ ~ 0.2/um). For this range of r, no experimental data 
are available. In contrast to the case of uniaxial systems, the motor behavior at the end 
of the filaments is crucial for the presence of jams in centered filament systems: If motors 
unbind actively at the filament ends (in the center of the aster), traffic jams are absent 
as shown in Fig. lEfa). In the regions with low motor densities, the density profiles for 
thermal (slow) and active (fast) detachment at the ends of the filaments agree perfectly. 

The main effect of the mutual exclusion is that density profiles get more and more 
flat with increasing motor concentration in the system, see Fig. Cfa). This means that 
the power law profile is found only for small overall motor concentrations. The average 
current in the system exhibits again a maximum at an optimal motor concentration. The 
maximum occurs at a motor concentration, where the bound motor density becomes nearly 
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constant and the power law behavior is hardly identified. 
Motors moving outwards 

For motors moving outwards in an aster we obtain profiles as shown in Fig. |7fb). For 
small numbers of motors (and not too close to the boundaries) the bound density follows 
the power law obtained from the linear equations. Now the exponent rj is positive but 
small. With increasing motor concentration, the profile of the bound density gets again 
more and more flat and the filaments become more and more crowded. As in the case of 
outward movements, however, the jams at the end of the filaments grow only very weakly 
and remain rather small (< 50£ ~ 0.5/im). The motor behavior at the filament ends is 
crucial for the presence of these jams also in this case, and jams are absent if motors 
unbind quickly at the filament ends, see Fig. IHfb). 

The new feature compared to the previous case is that the profile of the total motor 
concentration, which is rather flat for small motor concentration, develops a pronounced 
maximum in the center of the aster as the number of motors is increased beyond the 
optimal motor concentration, see Fig.Efb). This can be understood in the following way: 
If no ATP is added to the system, motors will accumulate in the center of the aster, 
simply because they bind strongly to the filaments, and, in the center, the number of 
binding sites per unit area is maximal. If ATP is added, motors are driven outwards by 
active directed motion. If now the number of motors in the system is increased, so that the 
motor movements are slowed down by the exclusion effect, the outward drift is suppressed 
and accumulation in the center is successively restored. 

Discussion 

We have presented theoretical results for the density profiles of molecular motors in arrays 
of cytoskeletal filaments. Motors were described as particles which move actively, i.e. in a 
directed manner, when they are bound to cytoskeletal filaments, but undergo non-directed 
diffusion upon unbinding from filaments. In addition, motor particles interact via mutual 
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exclusion. On the one hand, these models are designed to describe the generic behavior 
of the movements of molecular motors; on the other hand, the model parameters can be 
adapted to the transport properties of specific motor molecules. In general, these models 
involve certain simplifications compared to real systems. We have therefore tested a few 
modifications of the models to check that a more detailed modeling does not change our 
conclusions. 

Mutual exclusion of motors is obviously enhanced, if the motors carry large cargoes 
such as latex beads or vesicles. Furthermore, micr otubules consist of 12-14 protofilam ents 
which correspond to 12-14 parallel tracks (see also iNieuwenhuizen et al.l . .2002, 2004). For 
the uniaxial geometry, we have performed simulations of lattice models for which these 
two features have been incorporated. In these latter simulations, motors which occupy a 
cubic volume of M 3 lattice sites move on a microtubule consisting of 12 protofilaments 
arranged in a tubular geometry. In the simulations we chose M = 3 and M = 5. These 
cargoes can mimic small vesicles with diameters of some tens of nanometers. In vivo, the 
cargo diameters lie between a few nanometers for a single protein or RNA molecule and 
hundreds of nanometers for a large organelle. In addition, these model cargoes attach to 
the filament only with one of their surface sites, which represents the motor. Therefore, 
those cargoes which are not bound to the filament have an additional rotational degree 
of freedom. The resulting density profiles are averaged over M subsequent lattice sites, 
because the unrealistic cubic shape of the cargoes and the rigidity of their attachment 
to the motors leads to an artificial sublattice structure in the crowded region. We then 
obtain density profiles that resemble the ones discussed above, but the value of the bound 
density in the crowded region is smaller, because a smaller number of motors can block 
the filaments. In particular, a motor bound to one protofilament also blocks binding sites 
of the adjacent protofilaments because of the steric hindrance induced by its large cargo. 

Finally, let us relate our results to experiments. We have determined profiles of the 
motor density and motor currents in uniaxial and centered filament systems. On the one 
hand, these systems are directly accessible to experiments in biomimetic model systems 
in vitro. Density profiles as discussed here have so far only been measured for the case 
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of centered or aster-like systems (Nedelec et al 



2001|) (see our footnote at the beginning 



of the section on centered filament systems for a discussion of the dynamics in these 
experiments). The latter experiment shows the power law profile which is obtained from 
the theory for low motor densities. Higher motor densities and the corresponding jamming 
behavior have not been explored in this experiment, but could be studied in the same way 
by increasing the overall motor concentration. For the latter case, our theoretical study 
makes detailed predictions for the density and current profiles which could be checked in 
such an experiment. In addition, it would be quite interesting to construct other filament 
arrangements and compartment shapes and to study the corresponding motor transport 
experimentally. 

On the other hand, we can also compare our theoretical results about motor traffic in 
closed compartments with experimental studies on motor traffic in biological cells, where, 
however, additional phenomena such as the dynamics of the filaments, the regulation of 
the motor activity and the presence of other cellular structures play also important roles. 
Using fluorescence probes, several groups have measured the density profiles of molecular 
motors in vivo. One particular interesting system are kinesin motors in fungal hyphae. 
These hyphae are tubular compartments which contain uniaxial filament systems. In one 



exp erim ent, strong 



hyphae (S eiler et al 



ocalization of kinesin has been observed at the tip of these fungal 



2000). The comet-like density profiles of these motors localized 
at the tip correspond to the case of low motor density in our model. However, this 
localization is only found for kinesin mutants lacking a certain regulatory domain, i.e. 
for motors which m ove actively, but which are not regulated by cargo binding (see also 



Verhev et al 



1998). The underlying regulatory m echanism is the deactivation of th e 



motor via folding of its tail if no cargo is bound to it (|Cov et al 



1999 



Seiler et al 



2000). 



The deactivated motors do not exhibit active movement along filaments and can diffuse 
back over larger distances. Further regulatory mechanisms have mainly been discussed for 
the case of axon s where the question , whet her and how motors are transported back is 



most prominent (Go ldstein andY ang 



of motors at the axon terminal (Dahlstrom et al 



2000 ). The me chanisms include local degradation 



199ll ) and backward transport by motors 
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of opposite directionality (jHirokawa et al 



1990 



1991) 



Very recently, another fungal kinesin was also shown to localize at the tip of the 
hyphae and to exhibit these comet-like profiles. In this case, larger motor concentrations 
were induced by increasing the level of expression of the corresponding gene. This leads to 



density profiles w ith regions of hi g 



expression level (Konzack, 



2004; 



i motor density which increase in length with increasing 



Konzack et al 



200.4 ). 



According to our model, these 



density profiles should represent growing traffic jams. It would be highly desirable to 
repeat these experiments in vitro. 

In summary, we have discussed theoretically the stationary density and current pro- 
files of molecular motors in uniaxial and centered aster-like arrangements of cytoskeletal 
filaments. In particular, we have explored the effects of exclusion and jamming which 
can be addressed in these systems by varying the overall motor concentration. The two 
types of filament systems, which we studied, exhibit different density profiles and different 
jamming behavior. For small overall motor concentrations, the profiles are exponential in 
uniaxial systems, but algebraic in centered systems except for the crowded region close to 
the filament ends. Increasing the overall motor concentration, the jammed region grows in 
the uniaxial geometry, resulting in the coexistence of large regions of high and low density 
of bound motors, while the crowded region remains small in centered systems, where larger 
overall motor concentrations lead to a flattening of the profile if the motors move inwards, 
and to the build-up of a concentration maximum in the center of the aster if motors move 
outwards. In addition, the jamming in the uniaxial systems is rather insensitive to the 
motor behavior at the ends of the filaments, while the latter behavior is crucial for the 
presence of jams in centered systems. 

Both geometries studied here mimic arrangements of filaments in cells and are ac- 
cessible to in vitro experiments. The predictions for both geometries can thus be tested 
experimentally. Some density profiles have been de termined experimentally. These profiles 



correspond mainly to the case o f low motor density ( Nedele c et al 



only one recent experiment ((Konzack . 



2004; 



Konzack et al 



2001 



Seiler et al.1 



2000) 



2005) addresses higher mo- 



tor densities — and are in agreement with our theoretical description. 
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Appendix 



A Theoretical methods 
Monte Carlo simulations 

We performed Monte Carlo simulations for the case of uniaxial arrangements of filaments, 
where filaments are located within a cylindrical tube-like compartments and aligned par- 
allel to the cylinder axis which we take to be the x-axis. We take the filaments to have 
the same length L as the tube, but we checked that shorter filaments lead to very similar 
results. The cylindrical tube with radius R is taken to consist of all 'channels', i.e. lines of 
lattice site parallel to the filament, with u = (y 2 + z 2 ) 1 / 2 < R and < x < L. Reflecting 
boundary conditions are implemented by rejecting all moves to lattice sites outside this 
range. Within the closed tube the number ./V of motors is fixed. Each Monte Carlo step, 
corresponding to a unit of the basic time scale r, consists of N Monte Carlo moves. At 
each move, a motor particle is chosen randomly and updated according to the random 
walk probabilities. 

Two-state model 

Our Monte Carlo simulations show that the stationary profiles of the unbound motor 
density depend only weakly on the coordinates perpendicular to the filaments. In order to 
determine the stationary state, we can therefore use a two-state approximation, in which 
all unbound 'channels' are treated as equivalent and the motors can be in two states, 
bound and unbound. The stationary state is then characterized (i) by the balance of 
bound and unbound currents, jb and j u b, respectively, as given by 



with < n < L and (ii) by the change of the bound current as a function of the spatial 
coordinate n arising from the binding and unbinding of motors which leads to 



j b (n) = <p(n)j uh (n) 



(A.2) 



3b{n) - jb{n - 1) = 7r a d/9ub(n) 1 - Pb(n) - ep h (n) 1 - p uh (n) . 



(A.3) 
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The latter equation expresses the fact that, in the stationary state, the sum of all outgoing 
current is equal to the sum of all incoming currents at any filament site n and corresponds 
to Kirchhoff 's first rule for electric circuits. Here, pb an d p u b are the local number densities 
of bound and unbound motors, respectively. (In the next step, we will express the currents 
j'b and j u b in terms of these densities.) The coordinate n along the filament is given by 
the spatial coordinate x along the cylinder axis and the radial coordinate r for uniaxial 
and radial arrangements of filaments, respectively. (ft(n) is a geometrical factor and will be 
explained below. The binding and unbinding rates have been rescaled in Ea. IA.31 e = 2e/3 
and vf ad = 2vr ad /3. 

In addition, we express the bound and unbound motor currents as functions of the 
motor densities. For the tube geometry, we use the convention that the bound motors 
move to the right (the case that they move to the left is then obtained via the reflection 
symmetry). The bound motor current is then given by 



where Vb is the velocity in the absence of other motors. In the presence of many motors, 
forward steps are only possible if the filament site in front of a motor is not occupied. The 
probability of a vacant site is given by [1 — pb] > which leads to the reduction of the current 
as a function of density expressed in Eq. IA.4I 

For radial arrangements of filaments, we have to distinguish inward and outward move- 
ments of bound motors. The bound motor current is given by 



j h (x) = v h p h {x)[l - p h (x + 1)] 



(A.4) 




(A.5) 



and 




VbPb(r + 1)[1 - Pb(r)] 



(A.6) 



for outward and inward movements 



respectively. 



In all cases, the diffusive current of unbound motors is given by 



3ub(n) = D uh p uh (n + 1) - p ub (n) 



(A.7) 
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with the diffusion coefficient L> u b of unbound motors. Note that the latter expression is a 
discrete version of the usual diffusive current D^dp^/dn. 

The geometrical factor (ft introduced in Eq. IA.2I describes the relative weight of the 
bound and unbound currents and is given by the number of filament 'channels' per non- 
filament 'channel'. In general, (ft is a function of the coordinate n. For Nf isopolar parallel 
filaments within a cylindrical tube, (ft is given by (ft ~ -kE? /Nf. In particular, for a single 
filament, (ft is given by the tube cross-section. Notice that, within the two-state model, the 
number of filaments appears only via this geometrical factor and leads to a rescaling of 
the accessible volume for the diffusion of unbound motors. In the case of centered filament 
systems, the volume available for the unbound diffusion depends on the radial coordinate 
r, and the geometrical factor (ft increases linearly with r. In this case, (ft is given by the 
ratio of the free surface (i.e. not covered by filaments) to the area covered by filament 
'channels' which leads to 

,, , 2Trrh-N f £ 2 2irrh , , . , 

*w = N y - w , - (A-8) 

where Nf is the number of filaments, I 2 is the cross-section of a single 'channel', and h is 
the height of the slab, into which the aster is confined. 

At the boundaries, x = and x = L, terms corresponding to currents through the 
tube walls have to be omitted in Eq. I A. 31 Together with the normalization condition 

L r i n 

\Pb(n) + <ft{n)p nh {n) = 

n=0 

which fixes the total number N of motors in the tube, these equations form a system of 
2L nonlinear equations for the 2L unknown densities pb(n) and p n b( n ) with < n < L. 
We have solved thi s system of non- lin ear equations numerically using Newton's method 



N 

[pb(n) + 0(n)p ub (n)] = — , (A.9) 



with backtracking ((Press et al 



1992). The advantage of the two-state approach over 



the Monte Carlo simulations is that it requires less computation time, so that larger 
systems are accessible. In addition, simulations take particularly long computation times, 
if unbound diffusion is fast compared to bound movement, D^/v^l ^> 1, which is the 
case for cytoskeletal motors without large cargoes. In this case, the basic time scale 
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r of the simulations is much smaller than the step time t s ~ r/(l — 7) of the bound 
movements because 7 is close to one. In contrast, in the two-state approach, the necessary 
computation time is independent of the parameter values. In addition, within the two- 
state approximation, the computation time for several filaments (arranged in parallel or 
in an aster) is the same as for a single filament. 

As we do not distinguish between the different non-filament 'channels' in the two-state 
model, we neglect depletion layers close to the filaments as we discuss in Appendix IU1 in 
some detail for the tube geometry. In addition, a mean field approximation is implicit in 
the relations for the bound motor current as given by Eas. IA.4MA. 61 However, a comparison 
of the stationary profiles from the two-state approach with simulation results obtained for 
the case of truly equivalent unbound channels, for which the two-state approximation is 
exact, shows very good a greement. We therefore conclu de that, in contrast to the open 



tube systems discussed by 



Klumpp and Lipowskvl ([200.11 ) , the mean field approximation is 



quite accurate for the closed systems discussed here. 

B Low density limit of the two-state model 

Some analytical results can be obtained for the non-jammed low density regions both 
in uniaxial and centered filament systems. For this purpose, we consider the continuum 
version of the two-state equations. 

Uniaxial systems 

The continuum two-state equations for uniaxial filament systems are obtained by expand- 
ing Eqs. El and up to second order in the lattice constant 3 which leads to 

*W1 - Ph) - D h ^ = Aa<^ (B.10) 

W bJ^A>(l - Pb) - D ^~Q^ = *adPub(l - Pb) - ePb(l - Pub) (B- 11 ) 



3 We expand Eq.|S|taken at positions x and x—1 and average the results to get an non-ambiguous result. 
This agrees with the result obtained by expanding the time-dependent equations. The expansion leads to 
Dh = vt,l/2, but within the continuum equations, we can also treat as an independent parameter. 
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with the boundary conditions jb = fbPb(l — Pb) — D^dp^/dx = at x = and x = L 
which express the fact that no motors enter or leave the tube. These boundary conditions 
imply, via Ea. IB.10l that also the unbound motor currents vanish at the boundaries. In the 
low-density limit, hard core repulsion or exclusion can be neglected. This is appropriate in 
the non-crowded region, where p\, <C 1. For simplicity, we also neglect the bound diffusion 
terms, i.e., we consider the case = 0. On the one hand, this can be understood as 
taking into account only the first non-vanishing terms in the derivation of the continuum 
equations. On the other hand, a comparison of numerical solutions of the continuum 
equations with and without this terms shows, that the precise value of is largely 
irrelevant for the solution, as long as the detachment rate is small, which, however, is the 
case for processive motors. In the low-density limit the equations become linear, 

VbPb = D uh ^^ (B.12) 

Vh ~dx = nadPuh ~ epb ' (B.13) 

and, in general, the solution is given by a sum of two exponential terms. One term, 
however, decreases exponentially with x and therefore contributes only close to the left 
boundary, where it ensures the boundary condition of vanishing current and leads to a 
larger initial slope of the density profile. For sufficiently large x, the solution is therefore 
increasing exponentially along the tube, 



Pb 



{x)^Ne x ^, (B.14) 



where J\f is a constant and 



e = - We ~ 9 (B.i5) 



1 + - 1 



The last approximation is valid for small v\, and is also obtained from our first approx- 
imation, Eq. |1] above, where we assumed that unbinding and rebinding are equilibrated. 
The unbound density is given by 

PM = ^Pb( X ) + ^ d -^=M(^ + ^-) e^, (B.16) 
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i.e., bound and unbound density are proportional in the low-density limit. The first 
term of the factor relating bound and unbound density is the one obtained in the case 
of equilibrated transitions between the bound and unbound states (e.g., from linearizing 
Eq. G3) , the second one is a correction which shows that binding and unbinding are also 
driven out of equilibrium if ^ 0. (Note, however, that this term is of order v?, since £ ~ 
l/i>bj so that up to linear order in radial equilibrium still holds.) This correction term 
is positive, thus the current of motors binding to the filament at a given site, 7f a dPub(a ; )) 
is larger than the current of unbinding motors at the same site, ep^(x), which is easy 
to understand, since the motors bound to the filament are driven away by the drift v^. 
For small driving velocity Ub, we can replace the local balance of binding and unbinding 
currents at a site x by the condition 

ep h (x + Vb/e) sa vrad/Oub^), (B.17) 

which states that motors binding to the filament at site x, move for a distance v^/e before 
they unbind at site x + v^/e. Inserting the solution given above, we can check that this is 
fulfilled for small v^/i. 

~ep h {x + v h /~e) = ee Vh/i ^p h (x) ~ (e + v h /£)p h {x) = * a dPub(x). (B.18) 

The fact that more motors attach to the filament than detach from it, indicates that 
this solution cannot be correct for all x. In a system without mutual exclusion, unbinding 
will be larger than binding to the filament only at the end of the filament. In that case, 
we can account for unbinding at the filament end by assuming that all motors that would 
have detached in the interval [L, L + v^/e] are forced by the boundary to wait at the last 
binding site of the filament until they detach. Therefore the density at the filament end, 
Pb{x = L), is given by 

1 f L +— p 
p h (x = L)~±J £ dxAfe x K = | Afe L ^ (W( e £ - l) « ^J\fe L ^, (B.19) 

where £ is again the size of the binding site and the last relation is valid for small velocity 
^b; for which the ansatz given in Eg. IB.lfl is justified. A comparison with simulations for 
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the case without mutual exclusion shows good agreement of the density at the last lattice 
site with the values obtained by this procedure. In reality, however, there is hard core 
exclusion and the present solution holds only as long as the bound density is sufficiently 
small and breaks down at a certain x because of the exponential increase of the bound 
density. 



Centered systems 

For centered filament systems, the continuum limit of the two-state equations for low 
motor densities leads to 

, n dp uh dp h 
VbPb = <PorD nh — h D b —- (B.20) 



dph n d 2 p h 
or or z 



(B.21) 



for both inward and o utward movements. I n the case D\> = 0, these equations are equiv- 
alent to those used by 



Nedelec et al 



(2001) to describe their experimental results. These 



equations lead to 



n , e dp h v h d 2 p h D h d 3 p h \ dp h 
v h p h = D uh (j) r — + —z- - —^r + D h -—. 

Wad or 7r ad Or 1 7r ad dr A I dr 



(B.22) 



To recover the asymptotic solution of 
terms of order r* 1 . We obtain 



Nedelec et al 



(2001J) we assume p^ ~ r v and neglect 



V = 7rad^b/(e^o-Dub) 



(B.23) 



which can be larger or smaller than zero depending on the sign of the velocity v\,. Note 
that the bound diffusion coefficient does not contribute to this asymptotic result. In- 
terestingly, neglecting terms of order r 11-1 is equivalent to the assumption that binding 
to and unbinding from the filament are balanced locally. Hence asymptotically, bound 
and unbound densities are related by p u b = (e/^adVb) i n contrast to the case of uniaxial 
systems, and p u b decays with the same power law as p^. 
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C Depletion layer 

In this appendix, we derive an analytical expression for the radial profile of the unbound 
motor density for the case of a single filament located along the symmetry axis of a 
cylindrical tube. We consider the linearized diffusion equations which are appropriate for 
the low density limit or the non-crowded region to the left of the 'traffic jam'. 
The balance of bound and unbound currents is given by 

dpb(x) f R d 

VbPb(x) = D h — — ^D uh 2-Kj^ duu—p uh (x,u) (C.24) 

and the unbound motor density fulfills the stationary diffusion equation with cylindrical 
symmetry 

/ d 2 d 2 1 d \ 
Aib 7^ + + -H Pub = (C.25) 

which holds for values of the radial coordinate u with £ < u < R with the filament 
radius Rf ~ £ and the tube radius R. The solution has to fulfill the boundary condition 
dpub/du = at u = R. The longitudinal boundary conditions are the same as in the 
two-state model. Binding to and unbinding from the filament are described by 

dpb n 9V , ~ vr^ 2 , « ln 0fi x 

V b-Q^r - Dh ~Q^2 = ~ e P h + ^ad— Pub(a;, u = £), (C.26) 

which represents the boundary condition for p u b at u = £. The separation ansatz 

p uh (x,u) = e^f(u) and p h (x) = Afe*K, (C.27) 

where M is a constant, leads to 

f( ) = ^fo Mu/QY^R/O-J^R/QYoiu/O 

J[ > ne 2 Mt/OYiiR/O - Ji(R/OYo(m) ' K ' ' 



where J$ and Yq are Bessel functions of order zero of the first and seco nd kind, respectively, 



and J i and Y\ are the corresponding Bessel functions of the first order (|Abramowitz and Stegun 



1984). For small u the radial profile behaves as f(u) ~ ln(it/£). 

The balance of bound and unbound currents, Eq. IC.241 yields the condition 



Vh = ^ + D uh %I(£/t,R/0 
4 tt£ 
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(C.29) 



with 



f(u) _ n _f R/ * dzz J (z)Y 1 (R/0 - Ji(R/OYo(z 



(C.30) 



from which the localization length £ is determined numerically. 



D List of symbols 

Z?b (one-dimensional) diffusion coefficient of bound motors 

Z) u b diffusion coefficient of unbound motors 

f(r) radial part of the concentration profile in a tube 

h height of the slab to which a filament aster is confined 

jb local current of bound motors 

jub local current of unbound motors 

Jb spatially averaged current of bound motors in a closed tube 

I lattice constant, given by the filament repeat distance 

L linear extension of the compartment, 

i.e. length of the tube or radius of the disk 

L* length of crowded region (traffic jam) 

n coordinate along the filament in the general case, 

n = x and n = r for uniaxial and radial filament arrangements, respectively 

N number of motors 

Nj number of filaments 

M normalization constant 

r radial coordinate in the aster geometry 

R radius of closed tube 

t time variable 

u radial coordinate in the tube geometry 

«b velocity of bound motor 

x spatial coordinate parallel to the filament 
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y, z spatial coordinates perpendicular to the filament 

a probability for a forward step of a bound motor per unit time r 

(3 probability for a backward step of a bound motor per unit time r 

7 dwell probability of a bound motor per unit time r 

e detachment parameter, e/6 is the detachment probability per non-filament 

neighbor site per unit time r 

e rescaled detachment probability, e = 2e/3 

r\ exponent of the asymptotic density profiles in asters 

£ localization length or decay length of the density profiles 

7r a d sticking probability for a motor hopping to the filament 

7f a d rescaled sticking probability, 7r a d = 2-7r a d/3 

p average motor density (bound and unbound motors) 

Pb density of motors bound to the filament 

p u b density of unbound motors 

r basic time unit, defined by r = t 2 /D u b 

t s step time 

(f) ratio of the number of filament channels to the number of non-filament 
channels; uniaxial arrangements: tube cross-section/filament number, 
asters: radius-dependent effective cross-section/filament number 

00 parameter for the ratio of bound to unbound channels in the case of asters 
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Figure legends 

Figure 1: 

(a) Molecular motors perform active directed movements characterized by the bound- 
state velocity vy, along a cytoskeletal filament. After unbinding from the filament, the 
motor undergoes non-directed Brownian motion with diffusion coefficient -D u b- As motors 
are strongly attracted to filaments, mutual exclusion of motors from binding sites leads 
to molecular traffic jams. We study stationary states for two geometries: (b) uniaxial 
arrangements of filaments in closed tube-like compartments and (c) radial or aster-like 
arrangements of filaments in closed disk- like compartments. 

Figure 2: 

Profiles of (a) the bound motor density and (b) the corresponding bound motor current 
jb as functions of the coordinate x along the filament in the closed tube for three differ- 
ent motor numbers N as obtained from Monte Carlo simulations. The tube has length 
L = 200£ and radius R = 25£. The transport parameters are a = 0.01 — 2e/3 ~ 9.93x 10~ 3 , 
P = 0, 7 = 0.99, e = 10" 4 , and vr ad = 1. 

Figure 3: 

Average current Jb of bound motors (filled circles) and traffic jam length L* (open circles) 
as functions of the total number N of motors in the tube. Geometry and parameters of 
motion are the same as in Fig. |2 The data points at N/L = 0.5,1.74, and 5 correspond 
to the profile shown in Fig. [2J 

Figure 4: 

Two-state model: (a) Profiles of the bound motor density pb as a function of the spatial 
coordinate x parallel to the filament and (b) average bound motor current Jb as a function 
of the number ./V of motors in the tube as obtained from the numerical solution of the 
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discrete two-state model. The chosen tube has length L = 1000 and radius R = 25. The 
parameters of motion are = 0.01, D u ^ = 1/6, e = 10 -4 , 7r a d = 1. The numbers of motors 
in (a) are (from right to left) N = 100, 200, 300, 400, 500, 1000, 1500, 2000, 2500, 3000, 3500, 
and 4000. 

Figure 5: 

Two-state model: Profiles of the bound motor density pb for active and thermal unbinding 
of motors at the filament end: A motor at the end of the filament detaches with the same 
probability as at any other filament site (thermal unbinding, solid lines) or by a forward 
step, i.e. with rate a (active unbinding, dashed lines). L = 600, R = 25, parameters of 
motion as in Fig. EJ and N = 200, 800, 1400, 2000, 2600 (from right to left). 

Figure 6: 

(a) Profiles of the bound motor density pb (thick line) and the radius-dependent unbound 
motor density p u b as functions of the spatial coordinate x parallel to the filament. The 
lines for the unbound density show the profile at different distances from the filament 
(u = 1 to u = 25, bottom to top in the part left of the jam and top to bottom right of 
the jam). The dashed line indicates the exponential ~ exp(x/37.4) as obtained from the 
linearized equations IC.24M(T.26I (b) Radial profile of the unbound density in the low den- 
sity region, x = — 40(o), -20(D), 0(o), 20(A), 40(x), and in the crowded region, x = 80(«). 
The solid line shows the analytical result given in Eq. IC.28I The simulation data at all 
positions in the low density region agree with each other and with the analytical result, 
which shows that the profile has the product form (jC.27|) . The transport parameters are 
as in Fig. [2j and the tube has length L = 201 and radius R = 25. 

Figure 7: 

Concentration profiles for motors moving (a) inwards and (b) outwards in aster-like arrays 
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of filaments as functions of the radia l coordinate r. Parameters are for motor complexes 



as described by 



Nedelec et al 



(2001J), see text. The numbers of motors are (from bottom 
to top) N = 10 4 (o), 10 5 (o), 10 6 (A), 10 7 (□), 10 8 (*), 10 9 (>). The profiles shown are 
profiles of the total motor concentration p ~ p^,/(j){r) + p n ^. Because of the logarithmic 



scale, discrete data points are only indicated for small r. 



Figure 8: 

Thermal versus active unbinding of motors at the filament end in centered systems: Profiles 
of the bound motor density pb as a function of the radial coordinate r for motors moving 
(a) inwards and (b) outwards in aster-like filament arrays. At the end of the filaments, 
motors detach with the same probability as at any other filament site (solid lines) or by a 
forward step, i.e. with rate a (dashed lines). The inset in (b) shows the region close to the 
filament ends where the profiles for the two cases differ. The parameters are as in Fig. 
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